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Role  of  Chemical  Effects  in  Daytime 
High  Latitude  Trough  Formation 


1.  INTRODUCTION 

The  phenomena  of  high  latitude,  F-layer  troughs  have  been  observed  and  studied  for  some 
time. 1,2,3  These  studies  have  dealt  with  both  daytime  and  nighttime  troughs.  Although  it  was  believed 
that  the  former  were  of  an  irregular  and  sporadic  nature,  recently  it  has  been  demonstrated4  that  the 
occurrence  of  daytime,  high  latitude,  F-reglon  troughs  are  not  an  occasional  phenomenon  but  a  rather 
stable,  prevalent  effect. 

As  far  as  mechanisms  responsible  for  formation  of  these  troughs  is  concerned,  there  is  fairly 
general  agreement  that  these  are  linked  in  some  manner  with  the  polar  plasma  convection  pattern.5 
More  specifically,  however,  it  is  still  not  clear  which  particular  mechanisms  contribute  to  which 
particular  troughs.  In  general,  several  mechanisms  have  been,  and  still  are,  suggested  as  possible 


(Received  for  publication  31  January  1989) 

1.  Muldrew,  D.B.  (1965)  F-layer  ionization  troughs  deduced  from  Alouette  data,  J.  Ceophys.  Res. 

70:2635. 

2.  Moffet,  R.J.  and  Quegan.  S.  (1983)  The  mid-latitude  trough  in  the  electron  concentration  of  the 

ionospheric  F-layer:  A  review  of  observation  and  modeling,  J.  Atmos.  Terr.  Phys.  45:315. 

3.  Roble,  R.G.  (1983)  Global  dynamic  models  of  the  earth’s  thermosphere  and  ionosphere,  ESA  J. 

7:405. 

4.  Whalen,  J.A.  (1987)  The  daytime  F-layer  trough  observed  on  a  macroscopic  scale,  J.  Geophys. 

Res.  92:2571. 

5.  Sojka,  J.J.,  Raitt,  W.J.,  and  Schunk,  R.W.  (1979)  Effect  of  displaced  geomagnetic  and  geographic 

poles  on  high-latitude  plasma  convection  and  ionospheric  depletions,  J.  Geophys.  Res. 
84:5943. 


candidates.6  The  two  most  notable  of  these  are  transport  per  se  and  a  chemical  effect  that  is  in  turn 
induced  by  transport.  In  the  former,  plasma  produced  at  one  location  is  transported  from  one  region 
to  another,  retaining  its  original  properties  to  some  extent;  in  the  latter,  the  convection  of  plasma 
relative  to  the  neutral  gases  results  in  changes  in  chemical  reaction  rates. 

In  this  report,  we  present  calculations  and  arguments  that  indicate  that  chemical  effects  are 
unlikely  to  contribute  substantially  to  the  formation  of  the  observed  high  latitude,  daytime  troughs. 


2.  ION  CHEMISTRY 


The  chemical  reaction  with  which  we  are  primarily  concerned  is 


0+  +  N2  NO+  +  N 


(1) 


with  the  reaction  rate  kj.  This  's  one  of  the  two  dominant  chemical  loss  reactions  in  the  region  of  the 
F2  peak,  where  trough  formation  is  observed.  The  other  major  reaction 

0+  +  02  02+  +  O  (2) 

with  reaction  rate  k2  =  lOkj  at  1000°K  has  a  neutral  reactant  density  n[02]  <0.1  n[N2],  More 
important,  the  reaction  rate  k2  changes  only  moderately  with  temperature  compared  to  Reaction  (1). 
(See  Figure  1,  based  on  the  expressions  developed  by  St.  Maurice  and  Torr.7  )  Thus,  we  will  confine  our 
attention  to  the  reaction  Eq.  (1).  as  the  most  likely  candidate  for  contribution  to  trough  formation. 

The  reaction  given  by  Eq.  (1)  does  not  decrease  the  overall  charged  particle  number  density. 
However.  NO+  +  e~  ->  N  +  O  has  a  rate  constant  of  the  order  of  2  x  10  7  cm3 /sec  so  that  the  time  constant 
for  dissociative  recombination  is  of  the  order  of  1  minute.  Since  this  is  much  shorter  than  the  time 
constant  for  the  reaction  of  Eq.  (1),  we  may  equate  the  latter  with  the  loss  of  charged  particles. 


3.  THE  EFFECTIVE  TEMPERATURE 


Reaction  rates  are  commonly  given  in  terms  of  an  effective  temperature1 


m, 


eff  ~ 


m 


i  +  Mn 


Tn+  m 


M, 


i  +Mn 


m(M  n 
T,  +  — 1 — 11 
1  m 


u“ 


i  +  Mn  3k 


(3) 


Here  m(  and  Mn  are  the  ion  and  neutral  reactant  masses  respectively,  T,  and  Tn  are  the  ion  and  neutral 
temperatures,  u  is  the  difference  in  ion  and  neutral  velocities,  and  k  is  the  Boltzmann  constant.  For 
reaction  ( 1)  mt  =  16  amu,  the  mass  of  0+,  and  Mn  =  28  amu,  the  mass  of  N2. 

A  variety  of  expressions  have  been  used  for  ^efT’  but  none  are  appropriate  to  the  problem  at  hand. 
The  reason  for  this  is  that  we  do  not  have  a  steady  state,  but  rather  a  dynamic  one:  through  corotation. 


6.  Grebowsky,  J.M.,  Taylor,  H.A.,  Jr.,  and  Lindsay,  J.M.  (1983)  Location  and  source  of  ionospheric 

high  latitude  troughs.  Planet.  Space  Set  31:99. 

7.  St.  Maurice,  J.P.  and  Torr,  D.J.  (1978)  Nonthermal  rate  coefficients  in  the  ionosphere:  The 

reactions  of  0+  with  N2.  02,  and  NO,  J.  Geophys.  Res.  83:969. 
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Figure  1.  Comparison  of  the  Temperature  Dependence  of  the  Two  Main  Reactions  Controlling  the 
Peak  0+  Density  in  the  F-Region.  The  calculations  are  based  on  the  formulas  given  in  St.  Maurice 
and  Torr.7 


a  given  locality  enters  the  polar  convection  pattern,  and  then  interacts  with  a  changing  ion  velocity. 
The  neutral  velocity  as  well  as  both  temperatures  are  then  determined  by  the  governing  differential 
equations.  Since  the  two  time  constants  involved,  that  for  change  in  moving  through  the  convection 
pattern,  and  that  for  energy  or  momentum  exchage  between  neutrals  and  ions  in  the  F-region  are  both 
of  the  same  order  (hours),  a  solution  of  the  differential  equations  is  required  to  find  the  effective 
temperature. 

We  imagine  the  ion  velocity,  Uj,  to  be  given  by  the  polar  convection  pattern  and,  in  accordance 
with  some  polar  convection  model8  focus  our  attention  on  the  center  of  the  pattern,  where  the  electric 
fields  and  velocities  rise  to  their  maximum  values.  This  portion  is  linear,  and  as  a  given  region 
corotates  under  the  pattern,  we  assume  the  convection  velocity  rises  linearly  at  a  rate  a  (m/sec2).  The 
three  equations  for  un  (neutral  velocity),  T,.  and  Tn  are  the  neutral  momentum,  ion  energy  and  neutral 
energy  equations: 

(4) 

un)2]  (5) 

un)2]  (6) 

In  addition,  the  expression  for  u,(t)  is 

u,(t)  =  at  (7) 


au 


dt 


f* 


11  =  "vnl  <un  "  ui) 


dTt 

miVln 

[3k  (Tn 

dt 

m,  +  mn 

dTn 

mnVni 

[3k(Tj 

dt 

m,  +  mn 

In  these  equations.  vln  and  vnl  are  the  ion-neutral  and  neutral-ion  momentum  transfer  collision 
frequencies,  respectively,  and  mn  is  the  mass  of  the  dominant  neutral  species  (atomic  oxygen  at  the 
F2  peak).  Note  that  here,  mn  *  Mn.  Eq.  (4)  is  the  momentum  equation  for  neutral  particles,  Eq.  (5)  is 
the  energy  equation  for  ions,  and  Eq.  (6)  is  the  energy  equation  for  neutrals.  Note  that  vta  and  vnl  are 
related  by 


nimiV|n  =  nnmnvnt 


(8) 


For  the  F-region  peak  and  topside,  the  relevant  collision  frequency  is  0+  (the  dominant  ion)  with 
O  (the  dominant  neutral),  and  the  dominant  collision  process  is  resonant  charge  exchange. 
Consequently,  the  collision  freqencies  vtn  and  vnl  have  a  temperature  dependence  approximately 
proportional  to  the  square  root  of  the  relative  temperature,  Tr  =  (T,  +  Tn)/2.  However,  we  have  found 
from  numerical  integration  of  Eqs.  (4)  -  (6)  that  the  effect  of  this  temperature  dependence  is  small  and 
can  be  neglected  for  present  purposes.  In  the  following  analysis  the  momentum  transfer  collision 
frequencies  are  assumed  to  be  constant. 

Define  v  =  v^  ,  x  s  2vt,  and  un(0)  =  0.  Solving  Eq.  (4)  using  Eq.  (7),  we  obtain 


8.  Sojka,  J.J.,  Rasmussen,  C.E.,  and  Schunk,  R.W.  (1986)  An  interplanetary  magnetic  field 
dependent  model  of  the  ionosphere  convection  field,  J.  Geophys.  Res.  91:11281. 
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“i  -  un  =  7  (l  "  e_X/2) 

Let  q  s  nii/rr^.  p  =  v^/v^  =  ^n^/n^  .  T0  =Tn(0).  K  h  m1a2/3kT0v2,  T,  =  T,/T0.  and 
Tn  =  Tn/T0  Then  using  Eqs.  (7)  and  (9)  in  Eqs.  (5)  and  (6),  we  obtain 

S1'  +  e  T)] 

tlf!  -  T„)  +  K  (1  -  2e->«  +  e'l] 

These  equations  are  exactly  soluble  in  terms  of  elementary  functions,  but  the  details  are 
cumbersome.  Here,  we  make  use  of  the  facts  that  g  =  1  and  p  »  1  in  the  F-region.  Then  taking 
Tn(0)  =  Tj(0)  =  1.  we  find 

Tjlt)  =  1  +  Kx  -  2K(1  -  e~T/2) 

T n(r)  =  1  +  Kx  -  4K(  1  -  e-T/2)  +  K(l-  e'T) 

To  lowest  non- vanishing  order  in  t, 

Tj(t)  =  1  +  J  x2 


Tn(x)  =  l+~ -t3  (15) 

The  difference  in  order  is  due  to  the  ratio  p  =  vln/vnl  being  so  large,  which  in  turn  is  due  to  the  fact 
that  the  gas  is  weakly  ionized. 

Putting  Eqs.  (9),  (12).  and  (13)  into  Eq.  (3)  and  setting  m,  =  mn  =  16  amu  and  Mn  =  28  amu,  we  find 


Teff(t(  1  +  33kT0v2 


[22vt  -  19  +  16e~vt  +  3e-2vt] 


The  quantity  in  brackets  is  monotonically  increasing  with  time,  so  the  maximum  value  of  Tefr  is 
attained  when  atM  -  VM,  the  maximum  value  of  velocity  in  the  convection  pattern. 

(Om  =  1  +33kToV2t~2  t22vt^  "  19  +  16e_VtM  +  3e'2vtMl 


Further,  with  regard  to  variations  in  v,  a  maximum  value  of  [teff]M  is  obtained  for  v  =  0: 
(~  \  14  m.V.2  nijVw 

(TeffjMM  =  l+33  -kVf  *  1+  04242  -FfT 


14 

33  kT0 


while  for  vtM  =  1, 


(~  1  miVM 

lTerfJM=  1  +0.2816-^ 
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To  convert  Eq.  (IS'  into  one  containing  the  electric  field,  we  assume  an  ambient  magnetic  field  of 
B  =  5  x  10"5  T  and  rl  .ain 

(Tefr)MM  =T0+  0.3265Em2  (Em  in  mV/m)  (20) 

Substitution  of  this  expression  into  the  expression  for  the  reaction  rate  k]  given  by  St.  Maurice 
and  Torr  gives  kj  as  a  function  of  EM.  This  functional  relationship  is  illustrated  in  Figure  2  for 
to=  1000  K. 


4.  CHANGES  IN  F2  PEAK  DENSITIES 


To  utilize  Eq.  (20)  we  need  values  of  electric  fields.  Although  values  of  fields  in  the  polar 
ionosphere  as  high  as,  and  occasionally  higher  than,  100  mV/m  have  been  observed,  such  values  are 
not  common.  Since  the  daytime  troughs  are  observed  over  a  broad  range  of  solar  activities,  we 
examine  models  that  give  typical  average  values.  In  this,  we  have  the  recent  work  of  Sojka  et  al8 
which  provides  a  Kp  dependent  field  model.  In  the  dawn  dusk  plane,  for  the  maximum  value  of  the 
electric  field,  EM.  they  find 

60  10  +  6. 5Kp 


em  = 


n  10  +  1 ,7Kp 


(21) 


where  the  units  of  the  field  are  mV/m.  The  value  of  EM  thus  ranges  between  ~20  mV/m  for  Kp  =  0  and 
-50  mV/m  for  Kp  =  7.  Thus,  for  the  largest  average  fields,  the  effective  temperature  for  T0  =  1000  K 
would  be  1800  K.  and  kt  would  be  -1.2  x  10~12  cm3/sec,  a  value  ~2  times  greater  than  the  "undisturbed" 
value  at  1000  K.  If  it  is  remembered  that  Eq.  (20)  is  the  maximum  value  to  be  expected,  then  factors  of 
-2  are  the  most  to  be  expected  in  average  enhancement  of  the  reaction  rate  of  interest.  As  will  be  seen, 
such  increases  lead  to  decreases  in  peak  F2  region  densities  much  less  than  those  actually  observed. 
Fields  of  100  mV/m.  on  the  other  hand,  will  produce  an  effective  temperature  of  -4200  K.  and  a  value  of 
kj  =  1 .6  x  10~‘ 1  cm3/sec.  which  is  -30  times  higher  than  the  value  at  1000  K. 

Changes  in  F2  peak  densities  induced  by  relative  ionic  motion  may  be  studied  in  several  ways. 
Schunk  et  al9  insert  the  ion  drift  velocities,  assumed  to  be  constant,  steady  state  values,  directly  into 
the  conservation  equations.  They  find  that  for  an  electric  field  strength  of  50  mV/m.  the  F2  peak 
density  is  reduced  by  a  factor  of  1.5  from  the  zero  field  value. 

Alternatively,  one  can  arbitrarily  vary  values  of  the  rate  constant  k,  and  compute  ion  density 
profiles  as  a  function  of  rate  constant.  The  values  of  kj  can  then,  via  Figure  1,  be  associated  with 
values  ofTcfr.  and  using  Eqs.  (19)  and  (20)  with  ion  drift  velocities  and  attendant  fields.  Such 
calculations  have  been  made  using  the  daytime  ionosphere  model  described  in  Decker  et  al10 *.  the 


9.  Schunk,  R.W.,  Raitt.  W.J.,  and  Banks,  P.M.  (1975)  Effects  of  electric  fields  on  the  daytime  high- 
latitude  E-  and  F-regions.  J.  Geophys.  Res.  80:3121. 

10.  Decker,  D.T..  Daniell,  R.E.,  Jr..  Jasperse,  J.R.,  and  Strickland,  D.J.  (1987)  Determination  of 

ionospheric  electron  density  profiles  from  satellite  UV  emission  measurements,  in  The  Effect 
of  the  Ionosphere  on  Communication.  Navigation,  and  Surveillance  Systems.  J.M.  Goodman, 

ed..  U.S.  Naval  Research  Laboratory.  Washington,  pp.  685-694. 
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Figure  2.  Rate  coefficient  for  0+  +  N2  -»  NO+  +  N  as  a  function  of  EM  (maximum  convection  electric 
field)  based  on  our  Eq.  (20)  and  the  formulas  of  St.  Maurice  and  Torr.7  The  calculations  assume 
Tn  =  1000  K. 


F-region  portion  of  which  is  based  on  the  model  of  Anderson. 1 1  Figure  3  shows  density  profiles  for 
three  values  of  kj,  5x  10~13.  2  x  10'12,  and  10~n  cm3/sec  corresponding  to  Tejr  values  of  1000,  2100, 
and  3600  K  or  electric  fields  of  0,  60,  and  90  mV/m,  respectively.  Noting  the  zero  and  60  mV/m  fields 
we  see  a  decrease  in  F2  peak  values  by  a  factor  of  -1.75,  in  rough  agreement  with  the  decrease  found  by 
Schunk  et  al.9  Even  using  the  reaction  rate  corresponding  to  90  mV/m  results  in  a  density  reduction  of 
only  a  factor  of  4. 

Whalen4  finds  decreases  in  the  daytime  trough  to  vary  by  factors  of  3  to  10.  Considering  that 
Eq.  (20)  represents  a  maximum  estimated  value  of  Tefy,  and  further,  that  even  using  this  maximum, 
model  calculations  do  not  produce  changes  of  3  to  10  in  peak  densities,  it  would  seem  that  the  chemical 
mechanism  can  make  only  a  minor  contribution  to  the  daytime  troughs  observed  by  Whalen. 


1 1.  Anderson,  D.  N.  (1973)  A  theoretical  study  of  the  ionospheric  F-region  equatorial  anomaly  1. 
Theory,  Planet.  Space  Set  21:409. 
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Figure  3.  Electron  Density  Profiles  for  Three  Values  of  the  Rate  Coefficient  kj  (0+  +  N2  ->  NO+  +  N).  The 
calculations  used  a  Jacchia12  model  atmosphere  with  an  exospheric  temperature  of  -1000  K,  and  the 
solar  EUV  fluxes  of  Heroux  et  al13  (for  F  10  7  =  120).  The  ionospheric  model  used  for  these  calculations 
is  described  in  Decker  et  al. 


12.  Jacchia,  L.G.  (1977)  Thermospheric  Temperature,  Density,  and  Composition:  New  Models. 

Smithsonian  Astrophysical  Observatory  Special  Report  375. 

13.  Heroux,  L.,  Cohen,  M.,  and  Higgins,  J.E.  (1974)  Electron  densities  between  1 10  and  300  km 

derived  from  solar  EUV  fluxes  of  August  23,  1972,  J.  Geophys.  Res.  79:5237. 
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